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We study the effects of quantum statistics on the counting statistics of ultracold heteronuclear 
molecules formed by Feshbach-assisted photoassociation [Phys. Rev. Lett. 93, 140405 (2004)]. Ex- 
ploiting the formal similarities with sum frequency generation and using quantum optics methods 
we consider the cases where the molecules are formed from atoms out of two Bose-Einstein conden- 
sates, out of a Bose-Einstein condensate and a gas of degenerate fermions, and out of two degenerate 
Fermi gases with and without superfluidity. Bosons are treated in a single mode approximation and 
fermions in a degenerate model. In these approximations we can numerically solve the master equa- 
tions describing the system's dynamics and thus we find the full counting statistics of the molecular 
modes. The full quantum dynamics calculations are complemented by mean field calculations and 
short time perturbative expansions. While the molecule production rates are very similar in all three 
cases at this level of approximation, differences show up in the counting statistics of the molecular 
fields. The intermediate field of closed-channel molecules is for short times second-order coherent if 
the molecules are formed from two Bose-Einstein condensates or a Bose-Fermi mixture. They show 
counting statistics similar to a thermal field if formed from two normal Fermi gases. The coherence 
properties of molecule formation in two superfluid Fermi gases are intermediate between the two 
previous cases. In all cases the final field of deeply-bound molecules is found to be twice as noisy 
as that of the intermediate state. This is a consequence of its coupling to the lossy optical cavity 
in our model, which acts as an input port for quantum noise, much like the situation in an optical 
beam splitter. 



I. INTRODUCTION 

The ability to create molecules from ultracold atoms by 
means of magnetic Feshbach resonances 0,0 or photoas- 
sociation Q has opened up exciting new areas of research 
from cold collisions physics to the study of strongly cor- 
related systems |4 H |g an d the BEC-BCS cross-over 
HL El EH • Both fermionic and bosonic atoms 
have been successfully converted into molecules whose 
condensation has also been achieved [30 03. More 
recently, Feshbach resonances between different atomic 
species [iall3 have been observed, leading to the poten- 
tial production of heteronuclear molecules. 

In this context the question arises how the quantum 
statistics of the composing atoms influence the molecule 
formation process. Both Wouters et al. 01 an d Dannen- 
berg et al. have studied the association of fermionic 
molecules from a boson-fcrmion mixture. They found 
that atom- molecule oscillations dominate the dynamics, 
similar to the case of two bosonic species. 

In this paper we analyze molecule formation in quan- 
tum gas mixtures of atoms of all possible combinations 
of statistics. We find that at the level of molecule num- 
bers the molecule formation dynamics looks very similar 
in all three cases and that differences due to the different 
quantum statistics only show up in higher order corre- 
lation functions. We are therefore led to calculate the 
full counting statistics of the molecular fields which is 
a computationally hard problem, but contains complete 
information about the molecule number and all its mo- 
ments. 

Measurements of higher order correlations are typically 
difficult so that, in spite of their value in characterizing 



matter wave fields, progress in that direction has been 
made only recently. Ottl et al. [2(J have determined the 
full counting statistics of an atom laser beam by single- 
atom detection within a hi gh- finesse cavity. Following a 
proposal by Altaian et al. |2ll |. noise correlation experi- 
ments in the spatial domain have been used to measure 
the second-order coherence function of pair-correlated 
atoms from molecular dissociation |22| and from an op- 
tical lattice in the Mott insulator regime |23j . 

We focus on molecule production in a two-step pro- 
cess consisting of both a ma gnet ic Feshbach resonance 
and photoassociation [24L l25l l26l| . The idea is to use a 
two-photon Raman process to transfer the closed-channel 
molecules - to which the free atoms are coupled through 
the hyperfine interaction - to a much more stable deeply- 
bound state, where they are not so prone to three-body 
losses. This method is similar to the recent experiment 
by Partridge et al. [2?]], in which optical spectroscopy 
was used to probe the two-body state of paired 6 Li atoms 
near a Feshbach resonance. One difference, however, is 
that we are considering here the Raman transfer inside 
an optical cavity that allows to select the final vibrational 
state of the molecules. 

We investigate the cases where the molecules are pro- 
duced from atomic condensates, from a mixture of bosons 
and fermions, as well as from fermions. In all cases the 
atomic sample is assumed to be at temperature T = 0, 
and in the case of fermions we consider both the case of a 
normal Fermi gas and of a superfluid system. We derive 
master equations that describe the molecule formation 
process in the presence of cavity losses. For small parti- 
cle numbers, these equations can be solved numerically 
in the single mode approximation for bosons and the de- 
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generate approximation for fermions. Thus, we find the 
full counting statistics of the molecules [2^| . 

We compare these numerical results to a perturbative 
analysis of the correlation functions and to mean-field 
calculations valid for short times. They show good agree- 
ment with the full numerical solutions of the master equa- 
tions in the limit of large particle numbers. 

We find that the qualitative features of the molecule 
production rates are similar in all cases. As has been pre- 
viously demonstrated e.g. for the case of four-wave mix- 
ing [23,|33, this is because the formation of superradiant 
collective states in fermions with Fermi energies Ep at 
times short compared to the dephasing time h/Ep leads 
to similar effects as the Bose enhancement for bosons. 

Significant differences arc, however, evident in the 
full counting statistics presented here. Atomic conden- 
sates and Bose-Fermi mixtures lead to coherent molecule 
statistics in the intermediate state, whereas normal Fermi 
gases lead to molecule statistics similar to those of a ther- 
mal field. The coherence properties of molecules formed 
from two superfluid Fermi gases is intermediate between 
the two. 

The second-order coherence of the deeply-bound mole- 
cules is particularly interesting, as this field is twice as 
noisy as that of the intermediate state. The additional 
noise finds its origin in the vacuum fluctuations of the 
cavity field, a situation somewhat analogous to that of 
an optical beam splitter, where vacuum noise is injected 
through the empty input port. 

The remainder of this paper is organized as follows. 
Section II discusses the case of molecule formation from 
condensed bosonic atoms and shows in particular that 
the optical cavity acts as an input port for quantum noise 
that increases the quantum fluctuations of the molecular 
field. Section III addresses the case of Feshbach-assistcd 
photoassociation from a Bosc-Fcrmi mixture that con- 
sists of a BEC and a degenerate Fermi gas and shows 
that under appropriate circumstances that system can 
be mapped to the previous situation of two heteronu- 
clear BECs. In section IV we consider the case of two 
degenerate Fermi gases with and without superfluidity. 
We conclude the paper with a discussion of our results in 
section V. 



II. TWO BOSE-EINSTEIN CONDENSATES 
A. Model 

We consider first the Feshbach-assistcd photoassocia- 
tion of condensates of bosonic atoms of species A and B 
into molecules. The atoms are prepared in a two com- 
ponent Bose-Einstein condensate. A sudden change of 
the magnetic field strength projects the scattering state 
of two free atoms onto the tunable molecular state. The 
ratio of open to closed channel admixture of this mole- 
cular state determines atom-molecule coupling and size 
of the molecular state. In entrance-channel dominated 




FIG. 1: Schematic of the atomic and molecular levels. The 
free atom pairs in the relative hyperfine state |0) can be taken 
from two BECs, a BEC and a Fermi gas or two Fermi gases, 
the latter with or without superfluid order parameter. Ac- 
cordingly, the molecules in state |1), |2) and |3) will be de- 
scribed by a single mode, if they are bosons, or by many 
degenerate modes, if they are fermions. 



Feshbach resonances a strong coupling between open and 
closed channel leads to molecular states which are pre- 
dominately in the open channel and thus have a substan- 
tial overlap with the wavefunctions of free atomic pairs. 
This explains the large coherent atom-molecule oscilla- 
tions such as in the Ramsey- type experiments 0] . In nar- 
row resonances outside the universal regime the molec- 
ular state is mostly in the closed channel which reduces 
its size compared to the very weakly-bound molecules in 
broad resonances but has some open channel admixture 
which ensures that the molecular state is not orthogonal 
to the free atom scattering states. For a detailed discus- 
sion of molecular production in Feshbach resonances we 
refer to 0, |H3 ■ In this paper we assume that the molec- 
ular state can be approximately described by a quantum 
field satisfying bosonic commutation relations. This pro- 
cedure is valid since the density no and the molecular 
size i.e. range of the molecular wavefunction ro satisfy 
n Q rl < 1. 

The free atom pairs with total spin configuration |0) 
are coupled to closed-channel molecules with spin state 
These closed-channel molecules are then coupled to a 
manifold of electronically excited molecular states \2 V ) by 
a classical light field of frequency ui and Rabi frequency 
f2„(R), where v labels the vibrational levels. These states 
are finally coupled via a cavity mode u(R) of frequency 
lu c to the electro-vibrational molecular ground state |3) 
with coupling strength g. Regarding the strength of the 
coupling between level |1) and |3) we restrict ourselves 
to the case k 3> g£l v / A„, where k is the lincwidth of the 
cavity and the detuning Aj, is defined below, i.e. we are 
not in the strong coupling limit of cavity quantum elec- 
trodynamics. The atomic and molecular levels together 
with the couplings between them are illustrated in Fig. 
□ 

At zero temperature and for sufficiently weak atom- 
atom and molecule-molecule interactions, excitations and 
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quantum depletion can be neglected. Here we limit 
ourselves to this case and some implications of atom- 
atom and molecule-molecule interactions are discussed 
below. Then the center-of-mass motion of atom pairs 
and molecules can be treated in a single-mode approx- 
imation with center-of-mass wave functions </>i(R) and 
corresponding center-of-mass energies e^. 

Denoting the internal energy of the molecular state \i) 
relative to that of the unbound atom pairs |0) by u>i, i = 
1, 2, 3, we assume that the detuning A„ = (u)2m — ^l) — 
lui w (uj2,v — ^3)— lu c s> \g\, so that we can safely adi- 
abatically eliminate the upper states \2„), leading to ef- 
fective two-photon Raman transitions between the states 
|1) and 1 3). The mean field shifts of each level have been 
absorbed in the energies w,; but in all that follows we will 
neglect the change in the mean field energies due to the 
changing populations of the various modes. In an inter- 
action picture, in which the state |3) has the energy of the 
two-photon detuning S = (lu 3 + e 3 ) — (uji + ei) — (uji — uj c ) 
and uj = lui + (ei — ca — es) equals the binding energy 
of the Feshbach molecule, Fcshbach-assisted photoasso- 
ciation is then described by the effective Hamiltonian 

IMIMEIIII 



in the master equation 



where 



H = 5ifilm 3 + H Q1 + ff 13 , 



H Q1 = a'rh\b A bBe wt + H.c. 



(1) 



(2) 



describes the hyperfine coupling of two bosonic atoms 
in state |0) to closed-channel molecules in state |1) with 
coupling strength a' . 
The second term 



His = x'rn\m 3 a + H.c 



(3) 



describes the Raman transfer of the molecules to a 
deeply-bound final state, together with the emission of a 
photon into the optical resonator, with coupling strength 
x' ■ In these equations, bA,bs are annihilation operators 
for atoms of species A and B, a is the annihilation oper- 
ator for photons in the single cavity mode, and rh\ and 
rh 3 are annihilation operators for molecules in the states 
|1) and 1 3), respectively. All the field operators satisfy 
bosonic commutation relations. In Eq. © the field that 
induces a virtual transition from the molecular state |1) 
to a virtually excited state |2) is treated classically. 

The dynamics of the cavity field, damped at a rate k, is 
described by the familiar master equation for a damped 
harmonic oscillator |39( . In the bad cavity limit, k 3> 
l^li Ix'IV^Vj with N the maximum number of molecules, 
it can be adiabatically eliminated 01 to give a master 
equation for the reduced density operator of the atom- 
molecule system, p = Tr cav j ty [w] . In this limit the two- 
photon detuning S does not influence the dynamics and 
we set it to zero for the rest of this paper. This results 



dp 
~dt 



a'rh 1 bAbBe lu ' t + H.c, p 



(4) 



+ 7 {mim\pm 3 rh\ — m 3 m\mirh\p + iJ.c.j 



where 7 = \x'\ 2 /k. The dissipativc part is of the Lindblad 
form idp/dt oc A* pA - Ap + H.C. with A = mirh\. 

The first term on the right hand side of Eq. (@J rep- 
resents the conversion of pairs of atoms in state |0) into 
molecules in state |1), and the second term leads to the 
amplification of molecules in state |3). Note that the 
rapid decay of photons from the cavity prevents their 
reabsorption in a |3) — > |1) transition, resulting in the ir- 
reversible transfer of molecules to a deeply-bound state. 



B. Molecule dynamics 

Typical experiments start out with all particles in the 
atomic condensates and no molecules, corresponding to 
an initial Fock state with Na atoms of species A and Nb 
atoms of species B. In the absence of losses the master 
equation |0J conserves the total particle number, 



dt 



(tia + n B + 2ni + 2n 3 ) = 0, 



(5) 



where ha and 115 are the numbers of atoms of species 
A and B, and n\ and 713 are the numbers of molecules 
in states |1) and |3), respectively. The evolution of the 
system can therefore be described on the basis 



\n,x) 



I n A/B 



N 



a/b ~ n,ni 



x, n 3 = x) , (6) 



where < n < minjA^jA^} is the total number of 
molecules and < x < n is the number of molecules in 
the final state |3). 

From the master equation Q we can derive equations 
of motion for the occupation numbers rii = (m). As 
usual for nonlinear systems, these lowest moments are 
coupled to a hierarchy of equations for moments of grow- 
ing order. Truncating this hierarchy by factorizing higher 
order correlation functions we find 



h 3 = 2(1 + n 3 )n 1: 
f>i = —2(1 + n 3 )ni 



-iae lLOT P ■ 



c.c. 



P 



-(1 + n 3 )P — iae %UJT \nx{nA + n B + 1) - nAn B ] 



fiB = iae lQT P + c.c., 



(7) 



where we have introduced the atom-molecule correlation 



P 



(m{b A bB) 



(8) 



a = a' /j, Cu = 07/7, and the time derivative is taken with 
respect to the dimensionless time r = 7t. 

Fig. |21 shows the result of numerically integrating the 
master equation |0J and the mean field equations Q by a 
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FIG. 2: Dynamics of molecule formation from two heteronu- 
clear BECs. Occupation numbers for the three modes ri3 
(solid line), ni (dotted line) and tia (dash-dotted line) are 
obtained (a) from the master equation Q and (b) from the 
set of mean field equations In both cases a = 10, Q = 0, 
and N A = N B = 40. 



standard fourth-order Rungc-Kutta method for a = 10, 
Co = and N A = N B = 40. 

The early stages of the dynamics are characterized by 
nonlinear oscillations between atomic pairs in state |0) 
and the closed-channel molecules in state |1) at a fre- 
quency of order ol^/NaNb- For longer times, though, 
the system is dominated by the two-photon Raman tran- 
sition from |1) to |3) and the damping associated with 
the optical cavity losses. 

While the details of the approximate results differ from 
the exact dynamics, there is good qualitative agreement 
between the two approaches, especially for short times. 
This is particularly true for N A ^ N B , in which case 
an excess of one bosonic species is always present. This 
leads to pronounced, more linear oscillations of the pop- 
ulation of atoms of the other species between the states 
|0) and |1), as predicted by the mean field equations Q 
and illustrated in Fig. [5] (b) . 

From simulations for varying particle numbers we infer 
that the dynamics obtained from the factorized moment 
equations converge to the full quantum dynamics. In 
the sense of this convergence Na = Nb = 40 is close to 
the "large particle limit", at least as far as qualitative 
features are concerned. 



C. Molecule number statistics 

An important feature of the effective Hamiltonian 
Eq. (|1I3J) is that we can solve the resulting master equa- 
tion (0} numerically and thus we are able to determine 
the molecule number statistics, that is, the probability 
P- (r) to measure j molecules in state \i) at dimension- 




FIG. 3: Number statistics for the chosed-channel molecules 
formed via a Feshbach resonance from a heteronuclear BEC. 
The parameters are a = 10, U) = 0, and Na = Nb = 40. 
Probabilities are cutoff at 0.3 for clarity. 



less time 



pU( T )=Tr[\n i =j)(n i = j\p(T)] 



(9) 



Fig-Hshows the time dependence of the molecule num- 
ber statistics for state |1), P- (t). For short enough 
times this distribution is very nearly Poissonian, and it 
remains so, until (ni) comes close to its maximum value, 
i.e. the molecular field can be considered coherent during 
this stage. 

To support this interpretation we calculate the second- 
order coherence function of the molecules in state |1) in 
the short time limit, 



,(2 



2) = lim g (2) (r) = lim 



<7ftt(0)lfoI(T)7fti(T)Thi(0)} 



(77i t 1 (0)mi(0))(m t 1 (r)mi(r))' 

(10) 

by expanding the expectation values in Eq. (|10f> around 
r = 0. We find 



7 (2) 



(b A b A b B b B b A b A b B b B ) _ (N A - 1)(N B - 1) 



{b A b B b A h 



N A N B 



(11) 

showing that the intermediate molecules in state |1) arc 
approximately second-order coherent in the limit of large 
particle numbers Na,Nb — ► oo. Moreover, we note that 
an excess of one species of bosonic atoms leads to the per- 
sistence of the initial second-order coherence over several 
oscillations between the states |0) and |1). 

Fig.0]shows the molecule number statistics in state |3). 
For short times, the associated second-order correlation 
function is found to be 



,(2) 



(b\b\blblb A b A b B b B ) 



{N A ~1)(N E 



i) 



{b\b B b A c 



N A N B 



;i2) 



5 



0.3 



0.2 - 




0.1 40 



the language of quantum optics [29J, |30j . 

The hyperfine coupling results, as before, in the for- 
mation of closed-channel molecules in state |1), which 
can then be transferred to a deeply-bound state via a 
two-photon Raman process that changes the momenta of 
the molecules by the difference of the momenta of the 
absorbed and the emitted photons, q. 

In these approximations the evolution of the system is 
governed by the effective Hamiltonian 

H = £]Tmt k TO3k + #oi+#i3, (13) 

k 

where 



FIG. 4: Number statistics for ground state molecules formed 
from a heteronuclear BEC. The parameters are a — 10, lo — 0, 
and Na = Nb = 40. Probabilities are cutoff at 0.3 for clarity. 



and 



H 01 =a'J2 mlJKbe^ + H.c. (14) 

k 



and exhibits the factor of 2 typical of a thermal field. The 
antinormally ordered correlation function of the cavity 
mode that appears in the first equality in Eq. Ijl2|l clearly 
shows that the loss of coherence suffered by the molecu- 
lar field during transfer from the intermediate state |1) to 
the deeply-bound state |3) finds its origin in the vacuum 
field of the optical cavity. The two-photon Raman tran- 
sition can be viewed as four-wave-mixing between two 
molecular fields and two optical fields where one optical 
mode is in its vacuum and thus the enhanced noise in the 
final molecular state is due to the spontaneous nature of 
the emission of the photon: The optical resonator acts as 
an input port that injects quantum noise into the molec- 
ular system, similarly to the shot-noise injected through 
the empty input port of an optical beam splitter. The 
enhanced noise in field |3) also manifests itself in the very 
broad number distribution in Fig. 0] 



III. BOSE-FERMI MIXTURE 
A. Model 

We now turn to the case of Feshbach-assistcd photoas- 
sociation from a Bose-Fermi mixture that consists of a 
BEC and a degenerate Fermi gas. We show that this 
system, in a degenerate model, can be mapped onto the 
case of heteronuclear BECs considered in the previous 
section. 

As before, the atomic condensate is treated as a 
zero-temperature condensate described by a single-mode 
bosonic field with zero momentum, background scatter- 
ing as well as three-body losses are neglected. The Fermi 
gas is assumed to be at T = 0. 

Differences in kinetic energy eventually lead to dephas- 
ing between the fermions. We focus on the short time be- 
havior (shorter than h/Ep) where the entire Fermi gas 
resonds cooperatively. This leads to a degenerate de- 
scription of the atoms, or homogeneously broadened in 



#13 = %' E ^lk^Sk+qO + H.c. 



(15) 



Here b is the annihilation operator for atoms in the 
bosonic condensate, /k is the fermionic annihilation oper- 
ator for atoms with momentum k, a the annihilation op- 
erator for cavity photons, and mik and 777,3k are the anni- 
hilation operators for fermionic molecules in the electro- 
vibrational states |1) and |3) with momentum k. The 
parameters 5, lu, a' and \' are the same as before. 
We proceed by introducing the pseudo-spin operators 



l lk/k> 
4k™lk - /k/k 
l lk" 2 3k+q, 
- rlk 777 lk -,7t k+q ^3k +q 



(*kf 



(16) 



as well as the total spin operators 



k • 



E- 

k 



s z 



rpZ 



k 

E' 



k- 



(17) 



in terms of which the effective Hamiltonians (|14|) and 
(fl"5|l become 



H13 



a'S+be lut + H.c. 
x'f+a + H.c. 



(18) 



Applying then twice the Schwinger mapping between 
an angular momentum operator and the creation and an- 
nihilation operators of two bosonic modes, 



S+ 
f+ 



m\b B , 



(19) 
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where &b, mi and 777,3 are bosonic annihilation operators, 
we can finally formally map the Bose- Fermi system under 
consideration onto the heteronuclear BEC of the previous 
section. Note, however, that the physical interpretation 
of the operators involved, and hence the physics, is dif- 
ferent in both cases. 



B. Molecule number statistics 

In the homogeneously broadened limit under consid- 
eration the atom-molecule coupling treats all momenta 
identically, so that the fermionic modes are always in a 
collective state, since they start in a totally symmetric 
state. 

Introducing the total occupation numbers of the states 
|1) and |3) as 



and 



"1 



"3 



J2 ™lk™lk 
k 



E^ 

k 



(20) 



(21) 



the evolution of the atom-molecule system can be de- 
scribed in the same basis as in Eq. © 

\n,x) = \n b / f = N b/ f -n,ni = n-x, n 3 = x) , (22) 

where, as before, < n < min{Ni,, Nf} and < x < n, 
N b and Nf being the initial number of bosons and 
fermions, respectively. It follows from the Schwingcr 
mapping that the matrix elements are the same for the 
Bose-Fermi mixture as for the heteronuclear BEC. 

Since all fermionic modes are treated on the same foot- 
ing, the molecule number statistics for a single fermionic 
mode of momentum k are simply given by 



P(n lk = !) = !- P(n ik = 0) 



W 
N f 



(23) 



and the molecule dynamics of a single fermionic mode 
|i)k are given by Fig. |3 with the vertical axes normalized 
to unity. 

The statistics of the total molecule number in states |1) 
and 1 3) are the same as their heteronuclear BEC coun- 
terparts and are shown in Fig. and Fig. For short 
enough times the excitation of the various momentum 
states of the closed-channel molecules are independent 
and the number statistics in state |1) is therefore 



P (i) = 
3 



N 



J 



[P(n lk = 1)] J [P(ni k = 0)] 



(24) 



a result valid as long as the molecule statistics for dif- 
ferent k are statistically independent. For large parti- 
cle number Nf, the binomial distribution converges to a 
Gaussian with mean (ni) and variance (ni) (1— (ni) /Nf), 
very similar to a coherent state, see Fig. |3J For longer 
times the various momentum states cease to be statisti- 
cally independent. 



IV. TWO DEGENERATE FERMI GASES 
A. Model 

We finally turn to the case of Fcshbach-assisted pho- 
toassociation from two normal Fermi gases of atoms of 
species A and B. Treating the Fermi gases in the homo- 
geneously broadened limit and the bosonic molecules in 
a single- mode approximation as in section [HJ the corre- 
sponding effective Hamiltonian is 



H 



8rh\ifi3 



01 



H 



13) 



with 



and 



H Q1 = a'm[ J2 fAMf B ,^e luJt + H.c. 



H13 = x'??4 7713a + H.c. 



(25) 



(26) 



(27) 



Here fA,k and /s,k are the fermionic annihilation oper- 
ators for atoms of species A and B with momentum k, 
a the bosonic annihilation operator for the single cavity 
mode, fiii and rh^ the bosonic annihilation operators for 
molecules in the states |1) and |3). The parameters 8, uj, 
a' and x' are the same as before. 

Introducing the pseudo-spin operators ^lj 



®k ~ (*k) ~ /l.k/i?,-k' 

K = ^ (/Jl,k/A,k + fb-kfB-k - 

and the total spin operators 



(28) 



(29) 



results in the master equation for the reduced density 
operator 



dp 
~dt 



a'm\S-e tuJt + H.c.,p 



(30) 



7 I mim\pm 3 rh\ 



rh3m\m\m\p + H.c 



B. Molecule dynamics 

Using the basis defined in Eqs. © and (|22|l . the action 
of the operators S~ on a state \n,x) is 

m\S-\n,x) = m\S-\S = N/2,m s = N/2-n) (31) 
= \/n - x + 1\/Jn + 1)(N - n)\n + l,x), 

where N = minjiVyi, Ng} is the maximum number of 
molecules, and Na and Nb are the initial numbers of 
atoms of species A and B. 



FIG. 5: Dynamics of molecule formation from a heteronuclear 
Fermi gas mixture. Occupation numbers for the three modes 
ri3 (solid line), ni (dotted line) and ua (dash-dotted line) are 
(a) obtained from the master equation (1301 and (b) from the 
set of mean field equations (1321 . In both cases a = 10, Q — 0, 
and N A = N B = 40. 

A semi-classical approximation can again be obtained 
by deriving equations of motion for the expectation val- 
ues of the occupation numbers and factorizing higher- 
order correlations. This gives 

n 3 = 2(1 + n 3 )m, 

h x = -2(1 + ra 3 )ni + (-iae iOT P + c.c.) , 

P = -(1 + n 3 )P + iae lOT Us + S~) + 2n 1 (5 z )) , 

^-(S+S-) = -2iae~ lQT P(S z - 1) + c.c, 
dr 

— Z ) =iae lOT P + c.c, (32) 
where we have introduced the atom-molecule correlation 

P=(m\S-). (33) 

Fig. [3] shows the evolution of the atom number and 
the populations of the two bosonic states as obtained 
from a numerical solution of the master equation (130(1 
and the approximate semiclassical equations of motion 
§2$ for a = 10, Q = 0, and N A = N B = 40. We again 
observe nonlinear atom-molecule oscillations between |0) 
and |1). The mean- field equations still describe the qual- 
itative features of the full dynamics, although the agree- 
ment is not as good as in the bosonic case |4^. Note 
that the slope of (tii(t)) near r = is different from 
the bosonic case above. This finds its explanation in the 
different molecule statstics to be discussed next. 



C. Molecule number statistics 

Fig. shows the molecule statistics for state |1). In 



FIG. 6: Number statistics for the closed-channel molecules 
formed via a Feshbach resonance from a heteronuclear Fermi 
gas mixture. The parameters are a = 10, u) = 0, and Na = 
Nb = 40. Probabilities are cutoff at 0.3 for clarity. 




FIG. 7: Number statistics for ground state molecules formed 
from a heteronuclear Fermi gas mixture. The parameters are 
a = 10, Q = 0, and Na = Nb = 40. Probabilities are cutoff 
at 0.3 for clarity. 



contrast to bosonic atoms, they are now reminiscent of a 
thermal state, a result corroborated by the second-order 
coherence function 

(2) = (S+S+S-S-) = 2(N-1) 
(S+S-) 2 N ' 

characteristic for a chaotic light field. This is because the 
fermionic atoms act as independent " radiators" , much 
like the atoms in a thermal light source, leading to the 
different slope of (jii(t)) near t = in the bosonic and 
fermionic case mentioned above. 

Fig. [7| shows the molecule number statistics in the 
deeply-bound state |3). As already discussed, the cou- 
pling of the molecular field to the lossy optical cavity 
provides an input port for quantum noise, so that the 
number fluctuations in that state are twice those in the 
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intermediate state Specifically, we find 



(aaa^a^) 



(s+s+s-s- 
{s+s-} 2 



2(N-1) 
2 " N ' ( 35 ) 



which leads for N — > oo and short times to g' 2 ' = 4. 

We conclude this section by considering two Fermi 
gases subject to an attractive interaction. The ground 
state of such a system is approximately given by the BCS 
state 



\BCS)=U(u k + v k s+)\0}, 



(36) 



with Bogoliubov amplitudes and and the atomic 
vacuum state |0). 

The second-order coherence g^ of closed-channel 
molecules in state |1) becomes in that case 



fJ 



(2) 



(A/vy 



((N) + (A/vn 



(37) 



where A = V J2k UkVk ^ s * ne § a P parameter, V is the 
attractive two-body potential between fermionic atoms, 
and (N) = J^k \ v k\ 2 the average particle number in the 
system. 

For state |3) we get 
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(2) 



2{A/V) 4 



((N) + (A/vn 



(38) 



which again differs from the result for state |f) by the 
cavity factor {aaa) a)) / ' (aa)) = 2. 

For A = we recover the results for a normal Fermi 
gas g (2) = 2 for state |1) and g (2) = 4 for state |3). 
As the number of Cooper pairs A/V becomes macro- 
scopic, which is the case for a superfiuid system, we re- 
cover the results for bosonic atoms g^ = 1 for state |1) 
and g^ = 2 for state |3). This demonstrates that super- 
fluidity results in a degree of coherence of the molecules 
intermediate between those of a BEC and of a normal 
Fermi gas. 



V. DISCUSSION 

In the single mode approximation for bosons and 
degeneracy approximations for fermions, the dynam- 
ics of Feshbach-assistcd photoassociation of ultracold 
molecules from atoms with different statistics can be 
treated in a unified framework. At the level of molecular 
populations, the molecule dynamics are to a large extent 
independent of the quantum statistics of the atoms. In 
the limit of homogeneously broadened fermionic systems, 
the descriptions of molecule formation from bosonic 
atoms and from Bose-Fermi mixtures can be mapped 
onto the same Hamiltonian with identical initial states. 
Although the interpretation of the results in both cases 
is different, with obvious fundamental differences at the 



detailed level of description of the molecular fields, their 
global properties, when summed over the momenta of the 
molecules, are therefore largely similar [TsUli^. 

The similarities in the generation of closed-channel 
molecules from bosons and fermions stem from the col- 
lective behavior of these systems, which are well known 
e.g. in the context of matter-wave four-wave mixing 
p3. l30l | . Bose enhancement is built into the state sym- 
metrization procedure for bosons, while the collective be- 
havior of fermions results from the constructive interfer- 
ences that occur when a collective state is built up by an 
interaction that cannot distinguish between the different 
particles. In addition, the two-photon Raman process 
that transfers molecules from the intermediate state |1) 
to their final state |3) does not distinguish whether the 
bosonic molecules are built out of bosons or fermions. 

Despite this analogy, though, the statistics of the re- 
sulting molecular fields exhibit a distinct signature of the 
type of atoms from which they are formed. This differ- 
ence stems from the fact that in case the molecules are 
formed from pairs of atoms in Fermi gases, these atoms 
act initially as independent " radiators" , much like in- 
dependently radiating atoms in a chaotic light source, 
whereas molecule formation from atomic Bose-Einstein 
condensates is a coherent process from the very begin- 
ning [28|. 

In the degenerate model the kinetic energies of 
fermions are neglected. They destroy the collective ef- 
fects that give rise to the similarities between bosons 
and fermions at times larger than h/Ep |29l. l3fjl|. Uys 
et al. ^3 compared molecule association in a degenerate 
to a non-degenerate model and found the short time be- 
havior indeed indistinguishable between the two, whereas 
the overall conversion efficiency is significantly reduced. 
Szymanska et al. |43j studied the short-time dynamics 
following an abrupt jump in the atomic interaction for 
the entrance channel dominated resonances observed in 
40 K and 6 Li. They point out that in these systems large- 
scale atom-molecule oscillations are precluded. We stress 
that, in contrast to our studies focus on Fcshbach 
resonances which are closed channel dominated. Nev- 
ertheless, for longer times, the dephasing between the 
fermions will destroy the collective effects and reduce the 
conversion efficiency. 

Atom-atom, atom-molecule and molecule-molecule col- 
lissions lead to several modifications of the theory as it 
has been presented here. Interactions between atoms 
in a BEC give rise to quantum depiction, that can be 
taken into account by means of the Bogoliubov approx- 
imation for the noncondensed part. Usually, these cor- 
rections are very small with the dominant physics well 
captured by the single mode approximation we use in 
this paper. For fermions repulsive interactions lead to 
very minor quantitative changes that are hard to distin- 
guish from excitations due to finite temperature. Attrac- 
tive interactions between fermionic species on the other 
hand give rise to the BCS phase transition and we have 
included them in the framework of the BCS mean field 
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theory. Atom-molecule interactions and interactions be- 
tween atoms of the two different atomic species may lead 
to phase seperation j^- For closed-channel dominated 
resonances however the background scattering is usually 
neglectable if the system is tuned close to the resonance 
[l8j . Inelastic collisions give rise to losses which can easily 
be incorporated in our formalism through a phenomcno- 
logical damping term, e.g. in Eq. Q. Calculations of the 
full counting statistics, however, become much harder nu- 
merically, as Eq. ceases to hold. It is worth noting 
that losses have no effect on our conclusions which are 
concerned with the short time properties of the molecule 
formation process. Thus, we are justified in omitting 
them. 

The second-order coherence properties of the final 
ground state molecules are further influenced by their 
coupling to a lossy optical cavity that acts as an input 
port for quantum noise, in a fashion reminiscent of the 
role of the empty input port in an optical beam split- 
ter. This process, which is dominated by vacuum fluctu- 
ations, results in an additional factor of two enhancement 
in the second order coherence , both for bosonic and 



fermionic gases. 

Future work will extend our model to take into ac- 
count non-condensed modes of the bosonic fields. This 
will allow us to study finite temperature effects as well 
as dynamical depletion of the condensate in the course 
of the molecule formation. Furthermore, we will use a 
more detailed description of the two-body physics. In 
doing so we hope to be able to study the effects of a 
time-varying detuning, which is necessary to adequately 
describe sweep experiments. 
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